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Abstract 

Entanglement swapping between photon pairs is a fundamental building block in schemes using 
quantum relays or quantum repeaters to overcome the range limits of long distance quantum key 
distribution. We develop a closed-form solution for the actual quantum states prepared by realistic 
entanglement swapping, which takes into account experimental deficiencies due to inefficient de- 
tectors, detector dark counts and multi-photon-pair contributions of parametric down conversion 
sources. We investigate how the entanglement present in the final state of the remaining modes 
is affected by the real-world imperfections. To test the predictions of our theory, comparison with 
previously published experimental entanglement swapping is provided. 

PACS numbers: 03.67.-a, 03.67.Bg, 03.67.Dd, 03.67.Hk, 42.50.Ex 



* the corresponding author: ascherer@ucalgary.ca 



I. INTRODUCTION 



Quantum cryptographic communication and quantum key distribution technologies have 
matured to a level sufficient for commercial applications. Yet, distance limits impact on 
their usefulness. To date, the only realistic proposals for long distance quantum cryptog- 
raphy are still based on optical systems. Light is an optimal candidate to be a carrier 
of quantum information because photonic quantum states are durable due to their gener- 
ally weak interaction with the environment and are conveniently manipulated by means 
of linear optics and photon-detection. Moreover, photons are the fastest and one of the 
simplest physical systems for encoding quantum information. The technological challenge 
is to establish transmission channels over long distances with a high signal-to-noise ratio 
using real-world optical fiber settings or free space. Long distance quantum communication 
(LDQC) is hampered by a significant loss of photons with distance traveled. In particular, 
loss of photons due to absorption during transmission in fibers is characterized by the expo- 
nential rule t — 10""'^/^° , where t is the transmission coefficient, d is the distance traveled 
and a is the loss coefficient of the transmission medium in units of dB. A further limitation 
to the development of quantum communication over long distances is a constant detector 
noise level dominating over the exponential decrease of the signal. In an effort to overcome 
these obstacles and range limits of LDQC, quantum repeaters [1, 2] or quantum relays [3-7] 
have been proposed, which comprise entanglement swapping [8] as a fundamental building 
block. 

In principle, quantum repeaters enable any distance to be achieved. The basic idea of a 
quantum repeater is to split the long distance quantum channel into shorter segments and 
to distribute entanglement between the end nodes of these segments. Then, after purify- 
ing the noisy entanglement for each segment, the entanglement is extended over adjacent 
segments by means of entanglement swapping. The purification procedure is repeated for 
the extended segments, and the whole protocol reiterated until high-purity entanglement is 
established between the end points of the channel. A quantum relay works in a similar way 
as the quantum repeater, but without the entanglement purification procedure and without 
quantum memories. This makes it much more feasible as compared to the repeater, but 
docs not allow achieving arbitrary distances [7]. With both schemes the signal-to- noise ratio 
can be appreciably increased. However, experimental realization [7] suffers from a number 
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of imperfections, including imperfect sources of entangled pairs and imperfect detectors. 

The impact of experimental deficiencies on the security of a quantum channel as well as 
on sifted and secret key rates in quantum key distribution (QKD) is of considerable rele- 
vance and has been the objective of a number of recent investigations. In [9] Brassard et al. 
showed that channel losses, a realistic detection process comprising detector inefficiencies 
and dark counts, and imperfections in the qubit source drastically impair the feasibility of 
QKD over long distances. In particular, the implications of using attenuated laser pulses 
instead of idealized single-photon on-demand sources were examined, and it was shown that 
unconditional security is very difficult to achieve in long distance QKD based on a BB84 pro- 
tocol [10] with such weak laser pulses. In the same work Brassard and coworkers obtained a 
more optimistic performance for QKD schemes based on parametric down- conversion (PDC) 
sources [11, 12]. The consequences of using probabilistic photon-pair sources (as realized by 
PDC) instead of (non-existing) single-pair on-demand sources for quantum communication 
including entanglement based QKD have recently further been investigated, see [13-15]. 

For LDQC employing quantum repeaters or relays it is particularly important to examine 
the issue of how the entangled quantum states after an entanglement swapping operation are 
affected by experimental imperfections. It is clear that due to these imperfections the actual 
quantum states deviate from desired Bell states and have to be described by some mixed 
states. The impact of transmission losses and detector inefficiencies as well as dark counts 
on the performance of quantum relays has been recently examined by Collins et. al. in [7]. 
However, the probabilistic nature of photon-pair sources has not been considered in their 
work. As will be explained in the next section, the probabilistic nature of PDC also involves 
the possibility of multi-pair generation. Depending on the "brightness" of the sources, the 
emission of two (or even more) independent pairs of entangled photons from the same PDC 
source at a time becomes a more or less significant event leading to faulty detection clicks 
and incorrect conclusions with regard to entanglement and correlations. Thus, the multi-pair 
nature of PDC sources impairs a high fidelity realization of entanglement swapping. The 
investigation of the issue as to what extent the inevitable multi-pair contributions of PDC 
sources impinge on the performance of quantum communication protocols based on entangle- 
ment swapping, is the main motivation for the research presented in this article. The effect 
of multi-excitation events in PDC in addition to detector imperfections and transmission 
losses on quantum repeater performance has been examined perturbatively in the context 
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of single-rail entanglement by Brask and S0rensen in [16]. Furthermore, the atom-light en- 
tangled states produced by Stokes scattering in the DLCZ-scheme [2] are very similar to 
the hght-hght entangled states produced in a non-degenerate PDC process, as both are two- 
mode squeezed states. In works by Jiang et al. [17] and Zhao et al. [18], which extend the 
original DLCZ scheme to dual-rail entanglement, the consequence of multi-excitation events 
has also been treated in a perturbative way. In our work we choose a different approach, 
which is non-perturbative and essentially simpler, as it uses only the very basic toolbox of 
quantum theory and the principle of Bayesian inference. 

In this article, we elaborate on how the entangled quantum states after entanglement 
swapping are affected by experimental imperfections, particularly the multi-pair contribu- 
tions of PDC sources. We provide a non-perturbative theory for realistic entanglement 
swapping with imperfect photon-pair sources as well as imperfect detectors. In particular, 
we incorporate the multi-pair nature of PDC, transmission losses, and detector inefficiencies 
and dark counts non-perturbatively. Our theory enables us to obtain a closed-form analytic 
solution for the resultant mixed entangled quantum states after a real-world (noisy) entan- 
glement swapping operation. To test our theory, we compare its predictions with actual 
experiments on entanglement swapping, that have been previously published elsewhere. For 
this purpose, we derive a closed-form expression for the probability of four-fold coincidences 
of four detectors, two for the Bell-state measurement and two for monitoring the remaining 
two entangled modes, one on each side, depending on variable polarization directions of an- 
alyzers. This result allows us to calculate numerically the visibihty of four-fold coincidences 
for arbitrary parameter values characterizing real-world imperfections of the PDC sources 
and detectors. Finally, we inspect how the entanglement present in the final state of the 
remaining modes is affected by the practical deficiencies. The analysis makes it possible to 
suggest the implications of the imperfections on schemes using entanglement swapping as a 
fundamental tool, and to optimize parameter settings for maximum performance. 

In addition to the imperfections of sources and detectors, further problems are encoun- 
tered in a realistic entanglement swapping process. Imperfections of temporal overlap of 
the light fields on a beamsplitter as well as spectral mode mismatch constitute inevitable 
practical difficulties deteriorating the performance. Moreover, the entangled photon pairs 
prepared by two sources are expected not to be of the same quality. All of these problems 
would complicate the analysis and are not taken into account in the present study, but are 
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planned to be included in our future work. Here, we would like to focus on implications of 
imperfect sources and imperfect detectors only. In this way, our considerations provide a 
very useful upper bound on the amount of entanglement after swapping. As for the influence 
of mode mismatch on the performance, we would like to refer the reader to models developed 
in [19-21]. 

Our work extends the previous challenge [9] to the security of QKD protocols using 
imperfect sources as well as the investigations in [7, 16-18]. Efforts toward establishing 
reliable transmission of quantum systems over arbitrary distances will benefit from our 
careful analysis presented here. 

This article is organized as follows. In section 11 we develop our theory of real-world 
entanglement swapping. A mathematical description of imperfect photon-pair sources and 
imperfect detectors is provided. Using these models we derive a closed-form solution for 
the quantum states after a realistic entanglement swapping operation. In section III we 
apply our theory to making predictions with regard to entanglement verification in terms of 
the visibility of four-fold coincidences and compare our numerical results with experimen- 
tal entanglement swapping. We proceed with a discussion of the impact of experimental 
deficiencies on the entanglement of the resultant quantum state. We conclude with a brief 
summary and suggestions for future research in section IV. 

II. A THEORY FOR PRACTICAL ENTANGLEMENT SWAPPING 
A. Physical situation and setting 

In this paper we develop our theory for entanglement swapping of photon-polarization 
qubits as an illustrative example. However, the main issues of the theory apply to any 
other photonic realization of qubits [22] and the results and implications are independent 
of the latter. The basic experimental situation is illustrated in Fig. 1. Two parametric 
down conversion sources emit photon pairs into spatial modes a, b, c and d, where a and b 
correspond to the first and c and d to the second PDC source. In the ideal-case scenario, 
one entangled photon-pair is emitted into the a and b modes and another one into the c 
and d modes. For entanglement swapping, a joined Bell-state measurement is performed 
on the b and c modes. This results in projecting the remaining modes a and d onto an 
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PDC PDC 



FIG. 1: Entanglement swapping of photon-polarization qubits, based on two imperfect parametric 
down conversion sources (PDC) and a Bcll-measuremcnt with four imperfect photon detectors. 
Four spatial modes are involved, labeled by a, b, c and d. Two modes, one from the first and 
one from the second source, b and c, respectively, are combined at a balanced beam-splitter (B). 
The exits of the latter, denoted by b' and d , respectively, are directed to polarizing beamsplitters 
(PBS) and then detected at four detectors: one for the H and one for the V polarization of each 
of the c' and b' modes. This set-up forms an interferometric Bell-state measurement [23, 24]. The 
four detectors are inefficient photon detectors subject to dark counts. Their readout is denoted by 
(qrst). Given this readout we are interested in the entangled quantum state of the remaining a and 
d modes depending on experimental parameters characterizing the deficiencies of the experiment. 

entangled state, depending on the measurement readout of the Bell-state measurement. As 
a consequence, the photons in the outgoing a and d modes emerge entangled despite never 
having had an interaction with one another [8]. The entanglement previously contained in 
the a and b and the c and d photon pairs is swapped to the a and d photon pair. 

As explained in Fig. 1, in the case of polarization qubits a Bell-state measurement consists 
in combining the b and c modes at a balanced beam-splitter, then directing its output modes 
b' and c' to polarizing beamsplitters (PBS) and finally detecting the four alternatives c'jj, c^, 
6y and b'^ at four detectors. The readout recorded by these detectors is denoted by (qrst). 
Since a polarizing beamsplitter transmits horizontal and reflects vertical polarization, the 
readout "g" refers to mode c'jj, the readout "r" to mode c^, the result "s" to mode b'^ and 
the readout "t" to mode b'^. The range of values that can be assumed by these readouts 
depends on detector type and is clarified below. 

The main task of this paper is to provide a model for an implementation of practical 
entanglement swapping. The goal is to model a realistic experiment with practical deficien- 
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cies. In a real-world scenario, the PDC sources are imperfect, creating not exactly one pair 
of entangled photons, but a superposition of alternatives that also includes the vacuum, 
independent pairs of photon-pairs, and higher pair-number contributions. This has been 
investigated before up to second order, see e.g. [14, 15]. Furthermore, the detectors used to 
perform the Bell-state measurement are never perfect. They are usually inefficient to some 
degree, meaning that they sometimes do not detect existing photons. Detection inefficien- 
cies are even further increased due to transmission losses between the source and detector. 
In theory, the latter can always be effectively taken into account by being included in the 
detector inefficiencies. Moreover, the detectors are also subject to dark counts, meaning that 
they may click and indicate a detection event even if there are no photons incident into the 
detector. In this paper wc will make a distinction between photon-counting detectors and 
detectors which cannot discriminate photon numbers. In each case, the recorded readout of 
a Bell-state measurement with inaccurate detectors will be denoted by (qrst). In the first 
case, the measurement results "gr", "r","s" and "f can indicate any photon number n eNq, 
whereas in the second case they are records of yes/no events, namely either "at least one 
photon" or "no photons". 

In what follows, wc develop the basic ingredients of our theory. We begin with a the- 
oretical description of imperfect photon-pair sources. We proceed by providing a detector 
model that takes into account arbitrary detector inefficiencies as well as dark counts. Using 
a Bayesian reasoning approach we finally derive the resultant quantum state p^'"''* of the 
remaining a and d modes in a realistic entanglement swapping experiment conditioned on 
the readout (qrst) of an inaccurate Bell-state measurement. 

B. Modeling imperfect photon- pair sources 

Ideally, a photon-pair source would create exactly one entangled photon pair on demand. 
Such sources do not exist yet. Realistic sources are probabilistic generating photon pairs at 
random instances within those time intervals allowed by the (pulsed) pump laser, and occa- 
sionally emitting two or even more photon pairs, although the probability for higher order 

contributions is usually kept small. While other approaches exist, see e.g. [2], parametric 
down conversion is the most common way to produce entangled photon pairs. 

In PDC, a crystal with an appreciably large x^^^ nonlinearity is pumped by a laser field. 
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Each of the pump photons can spontaneously decay into a pair of identical (degenerate 
PDC) or nonidentical photons (nondegenerate PDC). The rate of pair generation using 
PDC is proportional to the x^^^ nonhnearity, the strength of the classical pump field and 

the interaction time. As shown in [25], a PDC process can be described and mathematically 
represented by a one-parameter SU(1,1) transformation of the vacuum state: 

T(7)|vac) = exp (^i'yK^^ |vac) , 7 e M . (1) 

Here, Kx is one of the generators {Kx, Ky, Kz} of the SU(1,1) group defined by the commu- 
tator relations 



, ky^ = -ikz , \ky, i^,] = ikx , [k,, k^ = iky . (2) 



For instance, in the case of type-I nondegenerate PDC, in which a pair of photons is created 
in the same polarization, and which we will consider throughout this paper, this generator 
is given by the following two-boson realization: 

= ^ + «v^v) , (3) 

where and fey are the annihilation operators corresponding to vertical polarizations of 
two different spatial modes a and h. 

As we can see, the resultant generated quantum state is not just a pair of photons, but a 
superposition of photon number states which particularly also includes the vacuum, pairs- 
of-pairs, and even higher order contributions. For small values of 7, the quantum state (1) 
after a type-I nondegenerate PDC can be approximated as 

T(7)|vac) ^ |vac) + i-^kx\^a.c) = |vac) + ^|0110) , (4) 

where the Fock notation jij/cZ) represents a state with i^j^ k, I photons in the a^, Uy, by, 6h 
modes, respectively. Please become aware of the chosen order "HVVH" . This convention 
will turn out to be convenient with regard to the description of entanglement swapping, as it 
coincides with the order of labels of the corresponding detection events in the readout (qrst) 
of the Bell-state measurement, cf. Fig. 1. The role of the vacuum state in the superposition 
(4) is to allow for the particular feature that the generation of the desired photon pair 
occurs at random instances of time. To be more precise, a photon-pair emission happens to 
be random within those periods of time, during which a pump laser field propagates through 



the crystal. That is, when a pump pulse is sent, it will either lead to down conversion or 
not. There is a high probabihty that PDC will not take place at all. The strong vacuum 
component implies this. It is for sure, though, that there cannot be down-converted photon- 
pair creations during time intervals between two successive laser pulses, i.e. when there is 
no laser field propagating through the crystal. The randomness of photon-pair production 
can be decreased by using a crystal with a larger x^^) nonlinearity or stronger pump fields, 
but this also happens at the cost of increased probability of the emission of multi-pairs of 
photons, which is disadvantageous and to be avoided as far as possible. On the other hand, 
multi-pair-emission events can never be completely excluded. 

For entanglement swapping, two PDC sources are required. As introduced above, the 
two difi'erent spatial modes of the first and of the second PDC source are labeled by a 
and b, and by c and d, respectively. In addition, the photons emitted by each source can 
have two mutually exclusive polarizations. We label them by H and V, corresponding to 
horizontal and vertical polarizations. Moreover, any symmetrical superposition of horizontal 
and vertical polarizations is possible. They have to be regarded as quantum alternatives 
and taken into account coherently. In this paper, we assume a type-I nondegenerate PDC 
for both sources. Furthermore, wc elaborate our theory for polarization entanglement. The 
preparation of polarization-entangled photon pairs can be experimentally realized using 
a pair of identical crystals stacked together such that their axes are orthogonal to each 
other, whereas the pump laser is diagonally polarized. Such a combination of two crystals 
effectively creates a PDC source producing two-mode squeezed states of the form given by 
Eq. (1) in each of the two orthogonal polarizations, with generators as given by Eq. (3) and 
the equivalent form for the H polarization. The total quantum state prepared by two PDC 
sources of this kind is then mathematically represented as: 



We parameterize the generated quantum state by x = 7/2 G M. Since x is usually much 
smaller than one, the value x^ is the photon-pair production rate of the PDC source, some- 
times also referred to as its brightness. We will also call x the efficiency of the source. In this 
paper we assume the same efficiency for both PDC sources. Too small values of x lead to a 
strong vacuum contribution, so that most of the time the sources do not emit any photon 
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pairs. As the value of x increases the pollution from higher down-conversions becomes more 
and more important, see also [26]. 

For the purpose of doing quantum optical calculations, it is convenient to express the 

state (5) in a normal- ordered loim. This can be done as follows. Following [27, 28], given 
two independent bosonic modes a and h, wc may choose a different basis of generators as 
compared to the basis in (2) in order to obtain a two-mode bosonic representation of the 
su(l,l) Lie algebra: 



K, := a 



^b^ , k_ ■= ah , ko:=]^ (a) a + Ph + l) . 



(6) 



The new generator basis {kQ,kj^,k_} of the su(l,l) Lie algebra satisfies the following 
commutator relations: 



[X_,X+] = 2Xo, [Ko,K^]^±K^. 



(7) 



According to [27, 28], the following normal-order decomposition formula holds for exponen- 
tial functions of the generators of the su(l,l) Lie algebra: 



exp [Q;+i^+ -I- cto-f^o + Oi-KJ\ — exp [74+ir+] exp [ln(Ao)-fi'o] exp [A_i^_] , 



where Aq, A± are given by: 



A, 

e 



{a±/e) sinh^ 



cosh^- (ao/2^)sinh^ ' 
[cosh^- (aio/2^) sinh^]"^ 



(8) 

(9) 
(10) 

(11) 



Using this decomposition rule we can derive the following special case which we need for our 
purpose: 



exp 



o)y^ ah^ 



exp 
X exp 
X exp 



uj{x) (a^a + 6^6 + l) 



where we introduced the definitions: 



0(x) itanhx , 
^{x) — -In [cosh x] 



(12) 

(13) 
(14) 
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Ix) = exp [4a;(x)] exp 
X exp 
X exp 



Each of the four factors in Eq. (5) is of this form and can be decomposed in this way. 
By doing so and using the fact that creation and annihilation operators corresponding to 
different optical modes commute, we arrive at: 

uj{x) (a.'ldn + aldy + blJ)^ + Pyby + c|jCh + + 44 + 44 

4>{x) (^^H^H + dyby + Ch4 + Cv4) |vac) . 

It is easy to see that the last two exponential factors acting on the vacuum state leave the 
latter unchanged. We are thus left with: 



(15) 



Ix) = exp [4a;(x)] exp (f){x) (a^4 + 44 + 44 + 44) 



vac) 



(16) 



This is the normal-ordered form of the quantum state generated by the two PDC sources. 



C. Modeling imperfect detectors 

We now present our theory of detectors which we apply to practical entanglement swap- 
ping. We begin with the description of ideal photon-number discriminating detectors and 
then stepwise allow for practical deficiencies. As a first step we take into account detector 
inefficiencies disregarding dark counts. In a second step we provide a detector model which 
also includes dark counts. By this means wc get a theoretical description of photon-number 
discriminating detectors that arc inefficient and subject to dark counts. Eventually, we will 
also have to acquiesce to the fact that most of the detectors presently available in labora- 
tories cannot discriminate photon numbers, but instead effectively measure whether there 
are no photons or at least one photon in a mode. These kind of detectors are referred to as 
threshold detectors [29] . Finally, we explain our Bayesian inference approach which enables 
us to calculate the resultant quantum state after entanglement swapping in a realistic situ- 
ation with imperfect detectors given the knowledge of it in the hypothetical ideal-detector 
scenario. 



1. Unit-efficiency photon-number discriminating detectors with no dark counts 



Ideally, we would like to have a unit-efficiency, photon- counting detector that exhibits 
no dark counts. Such a detector never clicks when there is vacuum and always clicks when- 
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ever there are photons present in a certain mode, and the strength of the dick provides 
information about the number of photons. Following [29], we refer to such a detector as 
an ideal, photon-number discriminating detector. It is mathematically represented by a 
photon-counting projection- valued measure (PVM): 

{n„=|n)(n|, 71 = 0,1,2,3,...} (17) 

with respect to the Fock state basis , n e No} of a certain mode. 



2. Inefficient photon-number discriminating detectors with no dark counts 



In practice, however, detectors always have a non-unit efficiency, meaning that even if 
photons are incident into the detector it has a finite probability not to trigger a click event. 
Throughout the paper, wc denote the efficiency of a detector by r], with < rj < 1, where 
f] = 1 means a 100% efficiency. Following [30], we model a detector efficiency rj by preceding 
a perfect, unit-efficiency detector with a beamsplitter possessing the transmittance r}. This 
is illustrated in Fig. 2, if we replace there the thermal state in the second beamsplitter 
port by the vacuum state. If p^^g is the input quantum state of the signal mode and C/bs (rj) 
represents the unitary evolution corresponding to the beamsplitter transformation in the 
Schrodinger picture, then the probability to detect q photons, where g e No, is given by 



nqrTr^gf) 



Ubs{v) (Psig ® I vac) (vac I) Ul^{r]) 



(18) 



where we first trace over the reflected mode then apply the PVM (17) and finally take the 
trace with respect to the transmitted mode incident upon the perfect detector. Hereby we 
made the agreement that "reflection" and "transmission" refer to the signal mode. 

For the purpose of this article, it is particularly important to expand on the special case 
in which the signal state is a photon number Fock state: p.ig = Given this input state, 

the conditional probability to measure q photons with a unit-efficiency detector would be 
p^=i(q\i) = Sqi, whereas the conditional probability with an ?7-efficiency detector {rj ^ 1) 
amounts to: 



Pr,{q\i) = Trt,, 



IlgTrrefl 



UBs{v){\^){A®\ysic){vsic\)UUr]) 



if i < g 



(19) 
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Thus, the conditional probabihty to detect q photons given that i photons are incident upon 
a non-unit efficiency detector is a BernouUi distribution, as intuitively expected. Clearly, in 
order to be able to detect q photons, the number i of incident photons must not be smaller 
than q, as we have not allowed for dark counts yet. 

Before we continue, wc would like to note that the detector efficiencies 7] of our theory 
are intended to take into account as much as possible all inefficiencies of the experiment. 
In particular, photon transmission losses, e.g., in fibers, filters, and other optical elements 
preceding a detector, can always be included in an effective efficiency of the detector. Thus, 
the efficiencies rj of our model are to be understood as effective efficiencies comprising the 
proper intrinsic detector efficiencies as well as all kind of other losses. 

3. Inefficient photon-number discriminating detectors with dark counts 

We proceed by including the possibility of dark counts. We simulate dark counts by 
assuming the environment to be not in the vacuum state but in a thermal state of the form: 



This density operator models a thermal source with an average photon number Tr(prn) = 
sinh^r and a pseudo temperature T — /ia;/(A;Bln[coth^r]). Thus, instead of assuming a 
vacuum state to be incident on the unused beamsplitter port, we combine our signal mode 
which is to be measured with a thermal state (20) at the beamsplitter with transmissivity 
f]. Our detector model is illustrated in Fig. 2. 

The physical motivation behind this model is the association that the origin of dark counts 
stems from stray photons incident onto the detector from the environment. Furthermore, 
since the imperfect detector to be modeled has a non-unit efficiency, not all of the photons of 
the signal mode nor all photons from the thermal radiation cause a click in the detector, but 
just a fraction of them. This is effective ly modeled by a beamsplitter whose transmittance 
Tj is intended to represent the non-unit efficiency, followed by a subsequent perfect detector 
with unit-efficiency and no dark counts. Let us stress, however, that the thermal radiation 
intended to be responsible for dark counts is entirely fictitious and does not have to corre- 
spond to a real, actually existing photonic field. Any source of noise causing dark counts 
of the detector, e.g. electrical noise inside the detector, etc., can be effectively simulated by 




(20) 
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FIG. 2: A model for an imperfect detector with efficiency r] and dark counts generated by a fictitious 
thermal background source. The signal mode and the thermal mode represented by the quantum 
states /3sig and /5t, respectively, meet at a beamsplitter with transmittance 77. One of its exits is 
directed to a perfect (ideal) detector (Did), the photons of the second exit port are discarded. 
The perfect detector is assumed to be a unit-efficiency photon detector with no dark counts. We 
will make a distinction of two cases. In the first instance we will assume the perfect detector 
to be photon-number discriminating. Later we consider the case where the perfect detector is a 
unit-efficiency threshold detector. 

an imaginary thermal photonic field coupling into the detector from the environment via a 
beam splitter. This has been proven for homodyne detection by Appel et al. in [31]. 

As will be shown in the appendix A, the probability for a dark count in a non- ideal 
threshold detector amounts to 



For any < 77 < 1, we can always find a pseudo-temperature T — T{r) to model any 
arbitrary value of p^^. The case r] — 1 seems to exclude a non- vanishing p^^ ^ 0. However, 
we may take the simultaneous limit 77 — > 1 and T — )■ 00 in such a way that any dark count 
probability is kept fixed. 

Using this detector model, we can calculate the conditional probability to detect q e 
No photons with a photon-number discriminating detector possessing arbitrary efficiency 77 
and any dark count probability pjc, given an input quantum state p^,^ of the signal mode, 
according to: 



(1 — 7]) tanh^ r 



(l_^)exp[-||] 
l_^exp[-||] 



(21) 



1 — 77 tanh^ r 




g|p,ig) = Trt,,„3 ngTr,,fl Ussir]) {Psis'^ PT)Uls{ri) 



(22) 



14 



Again, we use the convention that "reflection" and "transmission" refer to the signal mode. 
Furthermore, throughout the article, the subscripts 77 and p^c are used to express dependence 
on detector efficiency and dark count probability. For the purpose of this paper, it is sufficient 
to know the conditional probabilities for the particular input quantum states p^ig = 
i.e. Fock states. Wc arc therefore interested in p^p^^(g|i), with i,q & Nq, which is the 
conditional probability to detect q photons given that i photons are incident upon a photon- 
number discriminating detector with efficiency 77 and dark count probability pdc- Please note 
that now, due to dark counts, q may be greater than i. A derivation of these probabilities 
constitutes a significant technical part of the present work and can be found in the appendix 
B. The result reads: 

q 

(1 - vY G (i,q;v,&Ac) if ^ > (/ 

(23) 



t:tZi (1^) (l-^rG(.,g;,7,pJ iii>q 

1-?J(1-Pdc) 



ri'G{q,i;ri,p^,) ii q > i 



where 



b{rj, pdc) r^tanh r = 



1 + 



1 — 7] 

VPdc . 



(24) 



and the function G{-, •; rj, pdc) : Ng ^ M is given as follows. For A e No, > A, we define: 

1 2 



G{K,X;rj,pac) ■= Yl 



n=0 



K, — X + n 

K — X 



[b{V: Pdc)]' 



Fi ( — n, —A ; K — A + 1 ; 



V 



rj 



(25) 

and G{k, A; rj, pdc) '■— in the case k < X. Here and in what follows, 2-^i(-, ■);■;■) denotes 
the Hypergeometric function which is defined as 



n=l 



(7)n n\ 



(26) 



where (a)„ := r(a + n)/r{a) is the Pochhammer symbol, and r(-) is the Gamma function. 
Please note that ior q — i the two results in Eq. (23) coincide. 



4- Inefficient threshold detectors with dark counts 

If we go a step further and acquiesce to the fact that photon-number discriminating 
detectors are a technological challenge and their realization is still in its infancy (however, 
see [32, 33]), then, to provide a description of entanglement swapping of the most practical 
relevance, we have to consider threshold detectors [29]. Ideally, such detectors effectively 
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measure whether there are no photons or at least one photon in a mode. Following [29], 
we refer to a unit-efficiency threshold detector with no dark counts as an ideal threshold 
detector (ITD). Thus, an ITD is mathematically described by the PVM 

{no = |o)(o|, n>o = i-no} ■ (27) 

Inefficient threshold detectors with dark counts are contrived using the same detector model 
as above, but now with Did in Fig- 2 being an ITD instead of an ideal photon-counting 
detector. The relevant conditional probabilities are obtained in the same manner as above, 
using Eq. (22), but now with q being cither the event "no click" or the complementary event 
''click", corresponding to the PVM elements Ho or n>o, respectively. Again, for the purpose 
of entanglement swapping, we would like to know these conditional probabilities particularly 
for the signal input quantum states p^jg = The conditional probability of recording 

"no click" by a threshold detector with efficiency rj and dark count probability p^c given 
that i photons are incident upon it, can be directly obtained from the result (23) by setting 
q = 0. And the probability for the complementary event "click" is just one minus the latter 
probability: 

P'7,Pdc("nO Click"|i) = Pr,,pdc(9 = ON) 

= (1-Pac)[l-?7(1-Pdc)]' (28) 
Pv,Pac ( "click" |i) = 1 - Pr,,pdc ( "no chck" \i) 

= l-(l-pj[l-r/(l-pj]* (29) 

5. Bayesian updating based on evidence obtained by imperfect detectors 

In order to provide the resultant mixed quantum state of the remaining modes a and d 
depending on the result (qrst) of a Bell-state measurement on the c and b modes with imper- 
fect detectors including the presence of dark counts, we proceed in the fashion of Bayesian 
inference and reasoning. We first assume the notional ideal situation that the detectors used 
for the Bell-state measurement are unit-efficiency photon-number discriminating detectors 
with no dark counts. In this hypothetical case we know how to calculate the probability for 
a certain measurement readout {ijkl) of the perfectly accurate Bell-state measurement as 
well as the corresponding resultant pure quantum state, denoted by \^ijki), of the remaining 
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modes a and d after the measurement, using von Neumann's projection postulate [34]. We 
use this information as our hypothesis prior to observing evidence in the real experiment. 
The probability p{ijkl) is our prior probability of the hypothesis that the resultant quantum 
state of the remaining modes a and d after the Bell measurement is given by \^ijki)- 

Before we proceed, we would like to make the following agreement. Throughout the 
paper we agree upon using the letters q,r,s,t to denote the readouts of measurements using 
imperfect detectors, and the letters i,j,k,l to label results of hypothetical measurements 
employing perfect photon-number discriminating detectors. As for imperfect detectors, we 
differentiate between photon-number discriminating detectors and threshold detectors. In 
the first case q,r,s,t G Nq, in the latter case q,r,s,t ^ { "no click", "click"}. 

Given an actual detector readout (qrst) of an imperfect Bell measurement with inaccurate 
detectors including the presence of dark counts, we infer what an ideal four-tuple of detec- 
tors would have yielded, i.e., readout (ijkl), with probability -Pj^^f := p(ijkl\qrst). In the 
language of Bayesianism, after the evidence (qrst) has been observed by an imperfect Bell 
measurement, we update our knowledge with regard to the hypothesis according to B ayes' 
theorem, by means of which we can calculate the posterior probability of the hypothesis [ijkl) 
given the obtained evidence {qrst): 

_ p{qrst\ijkl)p{ijkl) ^^^^ 



Y.7,r,k',i'=aP((lfst\i'j'k'l')p{i'j'k'l') ■ 

Here, p{qrst\ijkl) indicates the conditional probability for obtaining the evidence (qrst) 
given that the hypothesis (ijkl) would have happened to be true, if the detectors of our Bell- 
state measurement had been ideal. It is important to realize that p{qrst\ijkl) is equivalent 
to the conditional probability of recording the readout (qrst) given that exactly i, j, k 
and I photons are incident onto the four non- ideal detectors of the Bell-state measurement. 
Hence, the resultant quantum state of the remaining modes a and d after recording the 
actual readout {qrst) at the inefficient detectors, is a mixed state of the form 

p''''-Y.pS\^^^'^i)i^^^'^i\- (31) 

ijkl 

In the next section we provide closed-form expressions for the conditional probabili- 
ties -Pjj^f depending on the experimental parameters x, rj and p^^. In the first instance 
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we assume imperfect photon-number discriminating detectors, i.e. detectors that reahze a 
photon-counting PVM but are in general inefficient (77 < 1), and have a non-zero dark-count 
probabihty p^c 7^ 0. As our next step we consider imperfect threshold detectors. Again, 
in the style of Bayesian reasoning we upgrade our probabilities of the density operator (31) 
using the evidence that the measurement outcomes g, r, s, t of a Bell-state measurement 
with threshold detectors can be either "no click" or "click". The resultant quantum state 
obtained in this way is the most relevant result, as it refers to the most common practical 
situation, namely entanglement swapping using ordinary inaccurate threshold detectors. 

D. Quantum state after entanglement swapping 

Using the tools presented in the previous sections wc are now in a position to derive the 
resultant quantum state (31) after a realistic entanglement swapping with imperfect sources 
and imperfect detectors. 

We start with the quantum state provided by the two PDC sources, Eq. (16). Suppose 
our four detectors of the Bell measurement on the modes b and c were perfect, i.e., had a 
unit-efficiency rj = 1 and no detector dark counts. Then, to give the ideal readout (ijkl), 
after applying the balanced beamsplitter transformation B = to modes b and c 

using the rule [35] 

(t)-K^:) ' (t)--(t)' 

and similarly for the vertical polarization, the resulting four- mode quantum state Ub\x) is 
to be projected onto the subspace corresponding to the projector 

The modes (^jj, Cy, by and b'jj are the output modes of the balanced beamsplitter, cf. Fig. 1. 
The operator Ub represents the unitary evolution corresponding to the beamsplitter trans- 
formation in the Schrodinger picture, which we use. Here and in what follows, \n) denotes 
an n-photon Fock state. It should be clear from the context to which mode it refers. Ac- 
cordingly, {\i){i\)c'^ represents a projection operator corresponding to the Fock state \i) of 
mode c'fj, etc.. The projection (33) followed by state normalization yields the following 
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post-measurement quantum state : 

\i'^'Hfvk''vf^)^\^,jki), (34) 



n?'l' Ub\x) 



n?'^ y y Ub\x) 

with the first factor being the Fock state with i, j, /c, / photons in the modes c'^, c'y, h'y, and 
h'jj, respectively, and 



\^i3kl) 




( V2)^+^+'=+'v/^!jIM! 5 S S ^^'^ W W U VA j ^^^^ 
x(at)-^(4)-(4)---^(4)^-— Ivac). 

Since the photons of the c and b modes are destroyed in the measurement process, we 
discard the first factor of Eq. (34). The second factor, \^ijki), is the resultant pure state 
of the remaining a and d modes: The corresponding probability of the hypothetical ideal 
measurement readout (ijkl) is given by: 

r, n -|2(i+j+fc+0 

p('m= i^l.,..,.pM = ^^^^^.^ ■ (36) 



H''-V'"V'"H 



In an actual experiment, however, detectors have an efficiency which is appreciably less 
than 100% and in addition exhibit dark counts. Given an actual detector readout (qrst) of 
an imperfect Bell measurement with faulty photon-number discriminating detectors charac- 
terized by efficiency rj < 1 and dark count probability p^^, we do not know the corresponding 
resultant post-measurement pure quantum state for the remaining a and d modes nor the 
probability of its occurrence. We can, however, calculate the posterior conditional proba- 
bility for any readout {ijkl) an ideal four-tuple of detectors would have yielded, i.e., the 
probability P^jkiiVipAc) = Pri,piS'>'jkl\qrst). As explained in the previous section, this is 
done using Bayes' theorem (30), according to which we achieve our goal if we know the 
conditional probabilities Pr,^p^^{qrst\ijkl) for all possible events {ijkl) e Nq. Since the four 
detectors are statistically independent from one another, these probabilities factorize into 
four terms: 

Pv,PaMrst\ijkl) =Pr,,pMi)Pv,pAr\j)Pv,pAs\k)Pri,pAi\l) , (37) 
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where each of these factors is given by the expression (23). The factorization of the condi- 
tional probabihties together with result (36) imply a factorization of the posterior probabil- 
ities of Eq. (30): 

Pt^fix, V, Pdc) = fHx, V, Pdc)/J(x, V, Pdc)/fe(x, V, Pdc)//(X, V, Pdc) , (38) 

where we now explicitely express their dependence on the experimental parameters x, rj and 
pdc- By inserting result (36) and our physical assumption (37) into (30), we have: 

Z^i'=oPv,Pdc\^r ) ^a^J^ \X) 
and equivalent formulae for /J (x, ?7, Pdc), fk{X:V:Pdc) and //(x, ?7, Pdc)- By inserting (23) 
into (39) we finally obtain the following exact closed-form solution: 

{ tanh^'(x)(l->?)'-'^G(i,g;r?,Wc) if i > O 

tanh2'(x)>?^('-^)(l-r,)^-'b^-'(r,,Pdc)G(g,»;r?,Pdc) jf « > ^ ^ ' 

with the common denominator 

9{q; X, V, Pdc) := tanh'^'(x) - vY-'b'-'iv, Pd.)G{q, t'; v, Pdc) 

i'=0 

oo 

+ tanh''{x){l-vy'-'G{t',q;r),p,^) (41) 

i'=q+l 

and similar expressions for f-{x,V,Pdc), fkix,V,Pdc) and //(x, ?7, Pdc)- Eqs. (38), (40) and 
(41) form our result for photon-number discriminating detectors. Knowing the probabilities 
PijkiiXi Vi pdc) for all possible hypothetical readouts (ijkl) with /c, / G Nq allows compu- 
tation of the mixed quantum state (31) for the remaining a and d modes after an imperfect 
Bell measurement with measurement result (qrst) using photon-number discriminating de- 
tectors with efficiency r) and dark count probability p^c- 

Let us consider the special case when dark counts are absent, i.e. pdc = 0. With Eq. (21) 
this leads to (1 — r^) tanh^ r — and therefore tanh^ r = as 77 7^ 1 in general. This in turn 
implies b{rj, pdc = 0) = r/tanh^r = as rj 0, and as a consequence also G{k, A; rj, p^c = 
0) = (^) . According to our general result (40) the absence of dark counts, pdc = 0, implies 
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i > q, again as b{r), pac = 0) must be zero, and we obtain: 



f!{X:V:P<ic = 0) 



Er=o r;")[(l-r?)tanh2x]' 
0[tanh^x]^-^(l-^)^- 



(42) 



(|l-[(l-»?)tanh2x]^ 



i.e., we get the result: 



f!ix, V, Pac = 0) = Q [(l-rj) tanh^ x] (l - [(1 - ^) tanh^ x] ) ■ (43) 
If we also let the efficiency rj go to unity, 77 ^ 1, we arrive at: 

lim fHx, V, Pdc = 0) = Sqi , (44) 

where 5qi is the Kronecker delta symbol. This is also what one should expect in the ideal 
case of perfect detectors with no dark counts. Please note, that with ^ 1 we have also 
implicitly assumed that there are no transmission losses between the sources and detectors. 

Acquiescing to the fact that photon-number discriminating detectors are still a tech- 
nological challenge, we now consider the more practical situation of threshold detectors, 
in which case the events (qrst) of the Bell measurement consist of readouts q, r,s,t e 
{"no click", "click"}. The conditional probabilities are given by (28) and (29). Calcula- 
tion of the functions //(x,?7, pdc) in Eq- (39) yields the following result: 

ii""°^"*'(X,^,Pdc) = [h{x,V:Pac)Y{l-HX:V:P,c)) , (45) 



I tanh^* X - (1 - pdc) [h{x, V, Pdc)f } 



i^"^'"" (X, V, Pdc) = i-p,, , (46) 

cosh X — 1 — T 

where we have introduced the definition: 

HX: V: Pdc) ■= [1 - 7^(1 - Pdc)] tanh^ x ■ (47) 

Eq. (38) together with (45 -47) form our result for threshold detectors. Please note, that in 
the ideal case, r} = I and p^c = 0, our result reduces to what we should expect: 

/r"""(x,^ = l,Pdc = 0) = 5,,o (48) 
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and 







if i = 



sinh"^ X 



(49) 



if i > 1 



oo 




^ytanh2^'x = l. 



(50) 



1=1 



As for the peculiar non-relevant marginal situation = and p^^ — 0, we note that in 

this case there can be no click events, so that f'^"^^"^" [j^^ ''7 = 0, p^c = 0) becomes a meaningless 
conditional probability. So it is not a flaw that the latter is not defined in this marginal 
situation. On the other hand, we get ''"'"'^ (x? ^ — 0, pdc = 0) = tanh^* %/ cosh^ x and 



Our derivation of the resultant quantum state (31) after an imperfect entanglement swap- 
ping is based on Bayesian reasoning and the physical assumption that the four detectors 
of the Bell-state measurement are statistically independent from one another implying the 
property of Eq. (37). The latter property involves a factorization of the probabilities Pfj^i 
into four terms corresponding to the four independent detectors, cf. Eq. (38). An obvious 
generalization of our results is to allow different efficiencies as well as different dark count 
probabilities for the four detectors. The derivation is very similar and the result reads: 

Pijkii X, {nA, {pdc,.}) = //(x, ?7i, Pdci)/J(x, ?72, Pdc2)/fe (X, ^3, Pdc3)//(X, ?74, Pdci) , (51) 

where rii, and pi^,, ^ = 1,2,3,4, denote the arbitrary and in general different efficiencies 
and dark count probabilities of the four detectors. The functions f^ix^Vi^^ Pdci/) are given 
either by the result (40) in the case of photon-number discriminating detectors or by (45,46) 
in the case of threshold detectors. 

III. COMPARISON WITH EXPERIMENTAL ENTANGLEMENT SWAPPING 

To cvahiatc the value of our model we test it against results of real entanglement swapping 
experiments [36-38]. In these experiments entanglement verification is accomplished either 
by observing the visibility of four-fold coincidences of four threshold detectors (two for 
the Bell-state measurement and two for monitoring the modes a and d, one on each side), 
obtained via variable polarization directions of analyzers, or by measuring certain correlation 
coefficients for polarization related to tests of the CHSH Bell inequahty [39] . Assuming the 



X/i=0 ft 



'no click' 



'"{x,v^ 0,pdc = 0) = 1. 
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resultant quantum state of the photons of the remaining modes a and d to be given by a 
Werner state, cf. [40], 

p^-* = F\^^)m + (U - , (52) 

where is some target Bell state, the visibihty is directly connected to the fidelity F ~ 
{^^Typqrst\^^T^ via the relation (cf. [38]): 

V"=(4F-l)/3. (53) 

Let us emphasize, however, that the assumption (52) is justified only for the post-selected 
quantum state plltLei, i-c, provided that click events are observed in both the a and d 
modes. This issue will be clarified below. The visibility is obtained according to V = 
(Max — Min) / (Max + Min), where "Max" and "Min" denote the maximal and minimal 
values of the four-fold coincidence rate as a function of some polarization angle, respectively. 
The relation between visibility and the correlation coefficient S of the CHSH Bell inequahty 
isS^2V2V, cf. [41]. 

To simulate a four-fold coincidence experiment and provide predictions for the visibility 

we proceed as follows. In accordance with the experimental situation of references [36, 38] we 
choose to consider four-fold coincidence events in which the result of the Bell measurement 
on the b and c modes corresponds — in the ideal case scenario — to a projection onto the 
Bell state 




According to the common argument on entanglement swapping, this would mean that the 
remaining modes a and d are left in the same Bell state, i.e. ]^~)a<i- To verify this entangled 
state by measurements on the a and d modes is the objective of the experiment of [36] . The 
experimental situations in references [37, 38] are very similar. 

Let us for a moment assume that our detectors and sources are ideal. Then, a projection 
onto the Bell state |^E'~)c6 is achieved whenever there are coincidence clicks of the two 
detectors for the c'^j and the by mode or, vice versa, of the two detectors for the c'y and 
b'fj modes. To understand this issue, we can exploit the fact that the Bell state |^~)c6 is 
antisymmetric under exchange of labels c and b implying fermionic statistics in the spatial 
behavior of the two photons, which means that they have to emerge from different ports of the 
balanced beam-sphtter (cf. Fig. 1). The other three Bell states are symmetric with respect to 
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exchange of labels c and b involving bosonic statistics in the spatial behavior, meaning that 
the photons will emerge at the same output port of the beamsplitter. Hence, in the ideal 
case of perfect sources and detectors, observing coincidences between two detectors on both 
sides of the beamsplitter is an experimental evidence for a projection onto the state |^^)cb 
and thus also a preparation of the state |^~)ad for the a and d modes. Only the coincidences 
of "c'jij and by" or "cy and b'jj" clicks are possible, but not of two clicks corresponding to 
the same polarization, i.e., "c^ and b'jj" or "cy and b'y\ since the polarizations in the Bell 
state |^~)c6 are anti-correlated. 

Yet, in a real experiment scenario the PDC sources and detectors are not ideal. In 
particular, wc have to allow for the rare but non-negligible faulty events due to emission 
of multi-photon pairs as well as dark counts, and also take into account transmission losses 
and detector inefficiencies, so that the actual quantum state of the remaining a and d modes 
after the suggested Bell measurement will deviate from \^~)ad- Within the setting of our 
model with imperfect threshold detectors, a non-ideal projection onto the Bell state \^~)cb 
is achieved whenever one of the following Bell measurement events (g, r, s, t) is obtained: 
(1, 0, 1, 0) or (0, 1, 0, 1). Here and in what follows, q = 1 means a click and q = that the 
detector does not click. 

Before we proceed we would like to comment on the following interesting observation. 
Even if the detectors of the Bell-state measurement were ideal (77 = 1, pjc = 0), they would 
never indicate a projection onto |^~)c6 and as a consequence prepare the Bell state \'if~)ad 
unless also the photon-pair sources were ideal. The multi-pair nature of the PDC sources 
precludes a projection onto the Bell state |^")c6 regardless of the quality of the detectors 
used for the Bell-state measurement. To understand this feature, let us consider the outcome 
(1010) of the Bell-state measurement on the b and c modes and assume ideal photon- number 
discriminating detectors. In this situation, according to Eq. (44), we have -Pj^^f = ^qi^rj^sk^ti 
and the quantum state (31) reduces to a single component, namely the pure state |$ioio)- 
According to (35) we obtain 



Thus, apart from the expected Bell state \ip~)ad, we get another term superposed to it, 
namely a superposition of two photons being in the a mode (in different polarizations of 
the latter) and no photons in the d mode and vice versa. This is also understandable 




(55) 
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intuitively. There are three quantum alternatives contributing to the event (qrst) — (1010) 
of the Bell measurement: (i) each PDC source emits exactly one photon-pair; (ii) the "first" 
PDC source produces vacuum and the "second" source two (independent) photon pairs; 
[iii) the "first" PDC source produces two photon pairs and the "second" source vacuum. 
The probability of each of these three alternatives is proportional to x"^, which also explains 
why the resultant normalized quantum state (55) of the remaining modes a and d does not 
depend on x- Only the alternative (i) leads to the desired Bell state |'0~)o(i, whereas the 
other two alternatives (ii) and (Hi) entail the second term in (55). For this feature see also 
[42]. 

Hence, entanglement swapping performed with PDC sources cannot herald Bell states 
in the outgoing modes with a 100% probability, even if the detectors used for the Bell 
measurement are perfectly ideal. If, however, we are interested in post-selected correlations 
of detection events of the a and d modes, only the Bell state |V'~)ad contributes to them. 
And as expected, the detection events will be perfectly anti-correlated. The second term in 
Eq. (55) does not contribute to four-fold coincidences and thus can be ignored in terms of 
post-selection. Furthermore, to fulfill the commonly used relation (53) between the fidelity 
and the experimentally observed visibility of four-fold coincidences, the former has to be 
calculated with the post-selected state. This means we have to project the quantum state 
-qrst Qj^^Q ^jjg subspace which corresponds to click events in both the a and d modes. Each 
of the outgoing modes has to have at least one photon. Introducing the projection operator 



■^-postsel 



>.,av - (|00)(00|)„^_„^) (lay,dH - mm)dy,dH) ' (^6) 

the post-selected quantum state is defined as: 

^qrst . n-postseiP npostsci (57) 

Tr [lIpoBtBeiP^''**IIpogtggi] 

The fidelity 



Fpostse. := (rKLir) . (58) 

fulfills the relation (53), (see [38]). 

Our scheme for entanglement verification via four-fold coincidence measurements is very 
similar to the experimental setup in [36]. It is illustrated in Fig. 3. We implement variable 
polarization measurements on both modes a and d by introducing polarization rotators 
(PR) into their spatial paths prior to polarizing beam-sphtters and threshold detectors. 
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The polarizations of the a and d modes are rotated by independently variable angles a and 
S. To avoid any confusion, let us agree upon the following meaning of the latter. The angles 
a and 5 stand for rotation angles of polarization vectors in the real space and not of Bloch 
vectors on the Bloch sphere. Neither do they mean rotation angles of the A/2-platcs, which 
are used in an actual experiment to cause polarization rotations. The relation between these 
three different meanings is as follows. A rotation by angle a of a polarization vector in real 
space corresponds to a rotation by angle a = 2q; of the Bloch vector on the Bloch sphere 
and to a rotation by angle ax/2 — a/2 of the A/2-plate. 

Polarization rotations prior to PBS separating horizontal and vertical polarizations, and 
photon-detections, are equivalent to polarization measurements in different bases in Hilbert 
space. The absolute angle of rotation in each of the modes determines the basis of the po- 
larization measurement. The polarization correlations of the entangled state should, in the 
ideal-case scenario, depend only on the relative angle between the two polarization rotators. 
We choose to rotate the polarization of mode a by a fixed angle a — n/A (i.e. -|-45°) and 
numerically calculate the probabilities for coincidences of detector clicks for measurements 
on the a and d modes, for varying angles S of polarization rotation of the d mode, given that 
the Bell measurement has yielded the result (1,0,1,0) or (0,1,0,1). We denote the four 
detectors for the measurements on the a and d modes by , D~ , and D^, correspond- 
ing to analyzing the a mode along the -|-45°-axis and — 45°-axis and the d mode along the 
variable polarization directions —S and +5, respectively. Furthermore, the four-tuple event 
(?2, ^2, S2, ^2) of the measurements on the a and d modes, cf. Fig. 3, corresponds to readouts 
of the four-tuple of detectors (D^, D~ , D'^ , -Dj"). 

In what follows, we provide the probability for four-fold coincidences using our model for 
practical entanglement swapping. Mathematically, the polarization rotators acting on the a 
and d modes are represented by the unitary operators 




(59) 



Ud{S) = exp i6Jd 



(60) 



with a, 5 e R and the generators of rotation 




(62) 



(61) 
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Here the angles a. and b are rotation angles of the Bloch vectors on the Bloch sphere and 
not of the polarization vectors in real space. As mentioned above, the former are related to 
the latter via a = 2q; , ^ = 25. 

Given an event {(i\r\S\t\) of an imperfect Bell measurement on modes h and c, the condi- 
tional probability that polarization measurements on the a and d modes would have yielded 
the result (i2j2^2^2), if the detectors had been ideal, is given by 



prob (i2j2A;2/2|giriSiti) = Tr 



X [U, 



■an -ay l^v jdn \ uaH aO-v l.dv ]dH\\ y 
*2 J2 ""2 ''2 /V2 J2 "'2 ''2 I ) ^ 



= E W^(«,^)/^--f(x,{77«},{p2l}) (63) 

iijikih 



PR a 




PDC 



PDC 



FIG. 3: A scheme for entanglement verification in an entanglement swapping experiment. A Bell 
measurement with imperfect threshold detectors is performed on the c and b modes. A read-out 
(qirisiti) is recorded. Polarization rotators (PR) parameterized by angles a and 6 are applied to 
the remaining a and d modes prior to polarization measurements using polarizing beamsplitters 
(PBS) and non-ideal threshold detectors. The readout of the second measurement is denoted by 
{(l2f2S2t2)- For a fixed polarization rotation of the a mode, a = 7r/4, the observation of certain 
Jour-fold coincidences for varying rotation angles < (5 < tt of the d mode reveals, via the visibility 
y, the strength of the polarization correlations between the a and d modes. 
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where we have used Eq. (31) and introduced the transition probabihties 

(64) 

The conditional probabihty to observe the event (g2'"2S2^2) with non-ideal, imperfect detec- 
tors, given an event (qiriSiti) of a non- ideal Bell measurement on modes b and c, is denoted 
and calculated as: 

QZZ% {x, WJ''}, WJ''}, {picl}, «, 5) := prob {q2r2S2t2\qins,h) 

= ^ prob(g2r2S2t2|«2i2A;2^2)prob(i2i2A;2^2kiriSi^i) 

- E ^P^pS(^2|^2)p,(2),,(.)(r2|J2)p,(2),,(^ X 

X ( E nS£{^^^)pS£{x,{v^h{^}) 

\iijikili 

00 00 

= E E (52K2) (^2b2) P^(2)_^(2) (.2|A;2) (^21/2) 

«2j2,fe2,'2=0 jljl,fel,il=0 

xW^; («> ( X, {r?!^^}, (65) 

For numerical simulations we need to know the transition probabilities (64). An 
explicit analytic expression is provided in the appendix C. This result together with 
Eqs. (51) and (45-47) enable us to calculate numerically the probabilities for four-fold co- 
incidences. As explained above, we condition on obtaining the readout (Ic'^, Oc'^, 1;,'^, Oy^) 
or (Oc^, Ic'^, Ofe/^, Ife/^) in the Bell measurement. Given either of these two events occurs, 
regardless which of them, we calculate the conditional probabilities of recording the events 
(1,0,1,0), (0,1,1,0), (0,1,0,1), and (1,0,0,1), using the four-tuple (L)+, D", D;^, £>+) of 
imperfect threshold detectors in the polarization measurement on the a and d modes, de- 
pending on the varying angle of rotation S. 

In order to compare our predictions with experimental entanglement swapping, we have 
to choose the same experimental parameter values {^t-^'*}; {v^^}: {pdcl}, {pfl} as in 
the experiment. It is important to emphasize that the detector efficiencies rj of our theory 
are to be considered as effective efficiencies that also include transmission and other losses, 
cf. Sec. lie 2. 

While we have developed our theory with entanglement swapping of polarization qubits 
as an intuitive example, it obviously also holds for any other realization of qubits, e.g. time- 
bin qubits [22]. In the following we compare the predictions of our theory with a recent 
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FIG. 4: Simulated probability of four-fold coincidence P4 depending on polarization rotation angle 
(5, for a fixed rotation angle a = 7r/4. Here a and S represent angles of polarization rotations 
of the a and d modes in real space and not on the Bloch sphere. The experimental parameter 
values are given by x ~ \/0706 « 0.24, r?^^^ = 0.045, 77^^] = 0.135, r?S5,3,4 = ^'^^ pic3,4 = 
1 X 10~^, p^^\ 2 = Pdci 2 3 4 = ^^ accordance with the conditions as found in the recent 

experiment [38] . The continuous sinusoidal curve refers to the probability of the click coincidences 
and D^" or and D^" , and the dashed sinusoidal curve displays the probability of the 
click coincidences "D~ and D^" or and D^'' , respectively, given that the Bell measurement 
yielded the result "iJ and V" or "F and if" for the c and b mode, respectively. The visibility is 

= 77.7%. 

experiment on time-bin entanglement swapping [38], which explicitly mentions all required 
parameters related to the PDC sources, transmission line and detectors. Fig. 4 and Fig. 5 
demonstrate results of our numerical simulations of four-fold coincidences with reference to 
experimental conditions of [38] . The conditions of this experiment are given by the following 
approximate values: x ~ a/0706 ^ 0.24, r][]l = 0.045, rji^} = 0.135, ^^2,3,4 = 0.04 and 
Pdc3,4 = 1 X 10~^, p^lJi^2 = Pd!i,2,3,4 = 3 X 10~^. As expected, we get two complementary sine 
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FIG. 5: Visibility V versus the square root of the photon-pair production rate X; for f?x 2 ~ 0.045, 
r/W = 0.135, 7?g 3^4 = 0.04 and p^^^A = 1 x 10-^ pi'\,2 = pi?i,2,3,4 = 3 x lO^^. These values 
refer to the experimental situation of [38]. The left image displays the functional dependence for 
very small x values, the right image reveals the behavior for appreciably large values of x- Foi' 
X ~ 0.24, which is the value found in the experiment [38], our numerical simulation (curve) yields 
the visibility V^^ = 77.7%. The experimental result V^^^ = (80 it 4)% of [38] is depicted by means of 
an ±4% error bar. A very good agreement is achieved. Our result also clearly demonstrates that 
the chosen photon-pair production rate belongs to a region with a rapidly decreasing visibility, 
thus being already far beyond its optimal value. 

curves that are 90° out of phase, one curve for anti-correlated polarization readouts ("-D^ 
and D'^" or and -D^") and another curve for correlated polarization readouts ("-D^ 
and D^" or and D^"), at the detectors for the a and d modes. As anticipated, the 
probability to detect anti-correlated polarizations for the a and d modes attains its maximum 
for 6 = a = 7T /4 and its minimum for 6 = a + tt/2 = 3tt/4. Complementary to this, the 
probability to detect correlated polarizations has its maximum for 6 = a + tt/2 = 3tt /4 and 
its minimum for 6 = a = 7c /4. This is the characteristic entanglement property of the Bell 
state |\&^)ad- The calculated visibility amounts to Kh = 77.7%. This is in respectably good 
agreement with the experimentally achieved visibility Kxp = (80 ±4)% of [38]. 

Fig. 5 reveals the behavior of the visibility as a function of the square root of the photon- 
pair production rate for detector efficiencies and dark count probabilities chosen as given 
above. It is interesting to observe that the photon-pair production rate of the PDC sources 
used in this experiment, x ~ 0.24, lies far beyond its optimal value. 
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FIG. 6: Visibility V versus the square root of the photon-pair production rate for various detector 
efficiencies r] and the fixed dark count probabihty p^^ = 1 x 10~^. The function ^,pdc(x) is plotted 
for six different efficiency values: rj = 0.025, 0.05, 0.1, 0.2, 0.3, 0.4, corresponding to the curves 
of lowest to highest visibility in all three diagrams. The left diagram displays the dependence of 
the function ^,pdc(/t) parameter r] for very small x values, the middle diagram shows the 

neighborhood of the maxima and the right diagram the behavior for high x values. The regions 
displayed in the first two plots are also included in the third diagram where they correspond to 
the steep increase of the visibility for very small x values. 

We have also compared the predictions of our model with the experimental data from 
the polarization entanglement swapping experiment reported in [37]. In this case, using x 
=0.05 (which can be calculated from reported count rates, transmission loss and detector 
quantum efficiencies), we find = 98%. Note that the experimentally obtained visibility, 
Kxp ~ 90%, is much smaller than the value obtained from our model. This can easily be 
explained by taking into account that detector noise and multi-pair emissions are not the 
only practical limitations that impact on the visibility. Other deficiencies include imperfect 
entanglement created by each source, even in the case where only one single pair is generated, 
and imperfect temporal overlap of modes b and c on the beamsplitter that realizes the Bell 
state measurement. We believe the latter to be the main limitation in this experiment. 
Hence, our model currently only provides an upper, yet useful bound on the visibility. 

To analyze how entanglement is affected by detector imperfections, we have calculated 
the visibility K),pdc(x) ^ function of x ^oi various detector efficiencies rj and dark count 
probabilities p^c- To keep the analysis simple, we have assumed the same efficiencies and dark 
count probabilities for all detectors involved in the experimental setting. The dependence 
on rj and pdc is illustrated in Fig. 6 and Fig. 7. It is interesting to observe that, according 
to our theory, for reasonable low dark count probabilities (p^c < 10~^) there is a region 
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FIG. 7: Visibility V versus the square root of the photon-pair production rate for various detector 
dark count probabihties and the fixed efficiency rj = 0.1. Plots of the function Vri^^^^ix) are shown 
for four different detector dark count probabilities: Pdc = 1 x 10^^, 1 x 10^^, 1 x 10^^, 1 x 10^^, 
corresponding to the curves of lowest to highest visibility in both diagrams. The left diagram 
displays the dependence of the function yr),pdSx) on the parameter p^c for small x values, the right 
diagram the behavior for high x values. 

of x-values such that the visibility can be made very close to 100% even for low detector 
efficiencies. Given any detector efficiencies and dark counts, our model makes it possible 
to provide the optimal photon-pair production rate in order to achieve high entanglement 
after entanglement swapping. This is a valuable result. An important conclusion that can 
be drawn from our investigations is the feature that high photon-pair production rates are 
counter-productive. One should not exceed values of about x = 5 x 10~^. Higher x values 
lead to a decreasing entanglement in the final quantum state after entanglement swapping 
and as such have an adverse impact. Yet, depending on the application, V may not be the 
appropriate figure of merit to optimize. 

First, in the here discussed proof-of-principle entanglement swapping experiments, the 
goal was to demonstrate a violation of the CHSH Bell inequality, rather than to maximize 
the violation. An important concern was thus to limit the time it took to take the data, 
generally many days. For the experiment [38] we referred to in Fig. 5, this has resulted in 
source brightnesses x exceeding the value required for maximum visibility. 

As a second example where the visibility is not the relevant quantity to optimize, let us 
briefiy consider QKD. Here, the figure of merit is the secret key rate, which scales in a non- 
trivial way with the visibility and with the number of detected coincidences [13]. For the 
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highest possible key rate, an optimal trade-off between the production rate of final entangled 
pairs, and visibility (i.e. amount of entanglement) has to be achieved. Hence, reducing the 
brightness of the PDC sources to a value that results in a maximum for V does in general 
not lead to a maximum secret key rate. 

IV. CONCLUSION 

Within the scope of our theory of practical entanglement swapping we have derived the 
actual quantum state of the remaining a and d modes depending on the result of a noisy 
Bell measurement with imperfect detectors and probabilistic sources. The main achievement 
consists in our ability to provide this quantum state for any given photon-pair production 
rate X, any given detector efficiency rj and any reasonable dark-count probability p^^. In 
deriving our results we have made a distinction between photon-number discriminating de- 
tectors and threshold detectors. The calculated quantum state allows us to make predictions 
with regard to any quantity of interest. In particular, we can calculate the entanglement 
of the remaining modes, parameterized in terms of the visibility obtained in coincidence 
measurements, depending on the experimental parameters V and p^c- 

Predictions of our theory have proved to be in close accord with experimental entangle- 
ment swapping. Our numerical simulations of certain four-fold coincidences between four 
threshold detectors demonstrate a remarkably good agreement with a recent entanglement 
swapping experiment reported in [38]. 

Furthermore, our theory provides a very useful functional relation between entanglement, 
quantified, e.g., by the visibility of four-fold coincidences, and the square root of the photon- 
pair production rate x- The latter is an essential experimental quantity in long-distance 
quantum key distribution based on protocols that employ quantum relays and repeater. 
The secret key rate, or equivalently the quantum bit error rate referred to as QBER, depend 
crucially on the photon-pair creation rate In our opinion, there sometimes is a tendency 
among scientists working on experimental QKD to aim at achieving brighter PDC sources 
due to the prevailed conception that higher photon-pair production rates lead to higher 
sifted-key rates. Here it should be realized, however, that it is not the source brightness 
itself that matters, but rather the production rate of final entangled photon-pairs after the 
swapping protocol, with their quantum correlations being as close as possible to that of Bell 
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states. Our results indicate that in some experiments the optimal value of x might eventually 
be exceeded. The investigations presented in section III suggest that, for a visibility of at 
least 97%, the photon-pair creation rates should not be chosen much greater than x — 0.05. 
Higher values decrease the entanglement present after swapping, due to undesired multi-pair 
events from the same source, thus implying an increasing QBER. Hence, our results allow 
us to suggest the implications of the imperfections on schemes using entanglement swapping 
as a fundamental tool. 

It should be emphasized, however, that, depending on the application, the amount of 
entanglement, as quantified by V, is not always the appropriate quantity to be made maxi- 
mal. In QKD, for instance, the figure of merit is given by the secret key rate R^^^, and the 
non-trivial dependence of the latter on x and the visibility V, cf. e.g. [13], suggests that the 
optimal X values with respect to QKD are not necessarily those which yield the highest vis- 
ibility of entanglement swapping. Nevertheless, we believe that the methods of the present 
work will prove very useful in finding the optimal photon-pair production rates with regard 
to achieving optimal secret- key rates in long distance, quantum repeater based QKD. As a 
first step into this direction we will have to generalize our results to several concatenated 
noisy entanglement swappings, and then examine the scaling properties with the number of 
these segments in the iteration. 

In our future considerations we also plan to investigate the issue whether the optimal 
X can be shifted to higher values by using photon-number discriminating detectors instead 
of threshold detectors for the Bell-state measurement. To understand the intuition, let us 
briefiy discuss the events where a coincidence detection in two threshold detectors in modes 
b' and c' is interpreted as a projection onto a Bell state. For threshold detectors, a frac- 
tion of these cases originates from two photons impinging on one detector, and one photon 
impinging on the other detector. These undesired but not identifiable events result in a 
reduction of the swapped entanglement, i.e. visibility. Unit efficiency photon-number dis- 
criminating detectors would allow identifying and discarding these events, which leads to a 
higher visibility for equal source brightness, or, conversely, allows increasing the brightness 
while keeping the visibility constant. And so would do imperfect photon-number discrim- 
inating detectors to some extent. This conjecture is worth analyzing in view of promising 
technological advancements in research on photon-counting detectors [32, 33]. Obviously, 
though, if the PDC sources are too bright, events where only one click is detected at each 
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output will be rare, and the production rate of final entangled pairs after the entanglement 
swapping operation will be low. 
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V. APPENDIX 



A. Dark count probability p^c 

Throughout the paper we characterize dark counts in terms of the experimental parameter 
pd^, which is defined to be the probability of a dark count event per detection window in a 
non-ideal threshold detector. Within our detector model, Eq. (21) provides p^c as a function 
of the parameter r , or equivalently the pseudo temperature of the thermal source (20). To 
derive this relation, let us consider our detector model as depicted in Fig. 2 and at first 
assume the ideal photon-detector Did to be photon-number discriminating. Furthermore, 
we also suppose the signal mode to be in the vacuum state: p^ig = |vac)(vac|. Then, the 
probability of detecting k dark count photons in the detector is obtained by applying 
the beam-splitter transformation C/bs(^) to the state |vac)(vac| 0Pt, followed by tracing over 
the reflected mode and projecting the transmitted mode onto the Fock state \k), and finally 
tracing over the latter mode. Hereby "reflection" and "transmission" refer to the signal 
mode. Thus, with the notation 

6(?7, r) := r/tanh^ r , (66) 

we have: 



prob( "k dark count photons") — Tr^^^^^ ^kT^vea t^Bs(^) (|vac)(vac| <S> Pt) ^bs(^) 

oo 



n=k ^ ' 



1—7] 

cosher V 77 
1 ^l-v 



k 00 



cosh r \ r] 



n=k 

^ r~ -ifc °° 
^%,r)J ^ 



v=0 



V + k 
k 



hM (67) 
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where we used fife = |^)(^| and the definition (20) of the thermal state px- Since < 
6(77, r) < 1, we may apply the Lemma 



E 

u=0 



u + k 
k 



[1 - qf+^ 



for < g < 1 , 



(68) 



which can be found in [43]. Hence, 



prob( "k dark count photons") — 



cosh r \ rj 



1-77 



1 - 6(77, r) 



1 - b{r],r) 



(69) 



To obtain the probability for a dark count event in a photon-detector that does not discrim- 
inate between photon numbers, we have to sum over all integers k ^ 0: 



Since 



we finally get 



prob( "k dark count photons") 
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(70) 
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cosh r \fc=o 
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(1 — 77) tanh^ r 



cosh^ 7- V 1 



(72) 



1 — 77 tanh^ r ' 

which is the result (21). Using this relation, we can express the function (66) in terms of 77 
and pdc rather than 77 and r: 



Kv,Pdc) := b{r),r{r),pac)) 



— 77 tanh^ r 



1 + 



1—7] 

VPdc , 



(73) 



which is Eq. (24). 
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B. Conditional probabilities P»7,pdc(9l*) 



The conditional probability Pr),p^^{q\i) to detect q photons given that i photons are inci- 
dent upon a photon-number discriminating detector with efficiency 77 and dark count prob- 
ability pd,. constitutes a major ingredient in our derivation of the quantum state after a 
realistic entanglement swapping operation. Here we derive its closed-form expression (23) 
using our detector model as depicted in Fig. 2. Assuming the signal mode to be in the Fock 
state Psig = according to Eq. (22) the conditional probability of interest is given as: 

Pr;,Pdc(9K) = Titrans HgTr^efl ^^^Bs (r?) ( K) | ® Pt ) t^Bs (r?) • (74) 

In what follows, we expand on the calculation of this expression. Using the relation (21), 
we expressed in (74) the dependence on the conditions of the thermal source in terms of the 
label pdc rather than the pseudo temperature of the source (or equivalently its parameter 
r). We will, however, initially use the notation pri,r{(l\'i) in the calculation below, and later 
show how the dark count probability p^c emerges in the final result. 

Let c, and e, denote the annihilation and creation operators of the signal mode and 
the thermal mode (fictitious environment), respectively. Initially, their common quantum 
state, denoted by ^3°g_T) is a product state, namely ^^"g^ = (S> Pt, where Pt is given by 
(20). We can express this input state as an operator functional of creation and annihilation 
operators: 

= [c, ct, e, et] = — \^{mer\y^c){yac\e-& . (75) 

cosh r ^!n! 

n=0 

Combining the signal with the thermal mode at a beamsplitter with transmittance r) leads 
to an entangled quantum state. In the Schrodinger picture, in which the photonic operators 
c, & , e, remain unchanged, the beam-splitter transformation is represented by the unitary 

operator Ubs{v) which acts on input quantum states of the beamsplitter. Its action corre- 
sponds to replacing the creation and annihilation operators in the functional g^^^ [c, c\ e, e^] 
by new operators according to the rule [35]: 



where 



B.-\ ;^^| (77) 
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is the beam-splitter transformation matrix consisting of the transmission and reflection co- 
efficients T — and R — \f\^^. Applying the beam-splitter transformation yields the 
following entangled quantum state: 

^i.T ^ uMQzAiv) 

= UM{\i){i\®PT)UL{v) 
tanh^"r 



cosh^ r 

n=0 



X |vac)(vac| — rjc+ \/rj€^ ^y^c — -^/T^-^e^ 

1 ^ tanh^" r ^ ^ ^ ^ ( \ ( i \ fn^ f n 
cosh^ r ^ i! n! 2^ 2^ 2^ 2^ L J U' J I J W 

n=0 i/=0 i^'=0 Ai=0 /i'=0 \ / \ / 



According to Eq. (74) the conditional probability Pr?,r(9K) obtained by ignoring the re- 
flected mode and projecting the transmitted mode onto the Fock state \q), and finally taking 
the trace of the resultant non-normalized state. This can be rewritten as: 

p^,,(g|^) = Tr[^:;^^^(flp®l-«)] . (79) 

By inserting (78) into (79) and realizing that in the Schrodinger picture, which we use here, 
the mode c in the equation (78) above corresponds to the transmitted and the mode e to 
the reflected mode, respectively, we get: 

1 X ^ tanh r \ ^ \ ^ x — ~ ^ ^ / t \ / t \ I n\ / n 



n=0 1^=0 i/'=0 /i=0 ;u'=0 ^ ' ^ ' ^ ' ^ 

X (-1)^+^ (^)2n+.+.'-,-,' (^^/T^) 
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—r-^j q\{i + n-qy. 



coshV (rr^) ' (1 - ")■ E [" t^""' r]" M,, n)]' .(80) 
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In the last step we have introduced the definition 



(81) 



It is important to recognize that fl{r),i,q,n) — unless i > q — n. This is intuitively 
understandable: the number q of detected photons in Di^ cannot be greater than the sum 
i + n of incident signal and thermal photons. To proceed, we have to distinguish two cases: 
i > q and i < q. It is easily seen that in the case i > q we have 

T) — 1 



n{rj,i,q,n) 
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(82) 



where 2-^i(-, •; S ") is the Hypergeometric function as defined in (26). In the second case, 
q > i, we obtain: 
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(83) 



We insert (82) into (80) and utihze the notation (66) to arrive at the following interim result 
for i > q: 
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For the q > i case, we insert (83) into (80) and obtain: 
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In the last step, after renaming the index of summation, ttt. — > ti, we made use of the 
permutation-symmetry property of Hypergeometric functions: 



2-^i(a, ^; 7; z) = 2-^1 (^, a; 7; ^) 



(86) 



In this way a partial symmetry between the two results (84) and (85) is achieved. Please 
note that the result (85) is valid also for g = i, as in this case the right-hand side expressions 
of (84) and (85) coincide. 

Wc would now like to replace the dependence on the parameter r of the thermal source 
by a dependence on the dark count probability. The latter emerges in the following manner. 
We replace the function b{7],r) by b{7], p^c) according to Eq. (73) and use Eq. (21) to make 
the replacement 

1 (1 -77)(1 - Pdc) ^g^^ 



cosh r 1 - 77(1 - pdc) 
Finally, by introducing the definition (25), we arrive at the result provided in Eq. (23). 



C. Transition probabilities Wj^^J^j'^^(a, d) 

In this appendix we provide an explicit analytic expression for the transition probabilities 
(64). Their definition involves unitary transformations representing polarization rotators for 
the a and d modes. According to Eqs. (59)-(62) the corresponding unitary operators are 
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given by: 



Ua{a) = exp 
Ud{S) = exp 



i 2" (^t^H + ayol) 



(88) 
(89) 



To find an expficit formula for W^^ltt 

{a, S) it is very useful to "disentangle" these expo- 
nentials using a generalized Baker-Cambell-Hausdorf formula. Let us define 



(90) 



These operators form a bosonic representation of the su(2) Lie algebra, as they obey the 
commutation relations of generators of the latter. According to [27, 28], the following 
decomposition formula holds for exponential functions of the generators of the su(2) Lie 
algebra: 
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We apply this decomposition formula to Eq. (88). In our case we have ^+ = ^_ = 1/2 and 
Oz = 0, yielding Li = 1/2 and 
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Thus, Eq. (88) becomes: 
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In a similar way we obtain a decomposition of Ud{S) in Eq. (89). The result reads: 
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Wc use these expressions in (64) to derive the transition probabilities. A straightforward 
but lengthy calculation yields the following result: 

with 

-I ii jl ki h / ■ \ / ■ \ / 1 \ / 1 

(^)n+.i+^.+^iV^T]7TTAJ^ 

'^^H+U+K+X,i2+j2 ^h+ji+ki+h ,i2+j2+k2+l2 



Y\ {k2 + k - jl - ki + ^ + n + ms) 



ms=l 



(100) 
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